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Abstract 

The dynamics of five dimensional Wilson line phases at finite temperature is studied 
in the one-loop approximation. We show that at temperatures of order T ~ l/L, 
where L is the length of the compact space, the gauge symmetry is always restored 
and the electroweak phase transition appears to be of first order. 

Particular attention is devoted to the study of a recently proposed five dimensional 
orbifold model (on S 1 /Z 2 ) where the Wilson line phase is identified with the Higgs field 
(gauge-Higgs unification) . Interestingly enough, an estimate of the leading higher- loop 
"daisy" (or "ring" ) diagram contributions to the effective potential in a simple five 
dimensional model, seems to suggest that the electroweak phase transition can be 
studied in perturbation theory even for Higgs masses above the current experimental 
limit of 114 GeV. The transition is still of first order for such values of the Higgs mass. 
If large localized gauge kinetic terms are present, the transition might be strong enough 
to give baryogenesis at the electroweak transition. 



1 Introduction 



Theories in extra dimensions are interesting alternatives to more standard four-dimensional 
(4D) extensions of the Standard Model (SM). In particular, there has recently been a lot of 
interest in theories with compact extra dimensions at the TeV scale [TJ . On the theoretical 
side they are, among other things, a possible framework to solve the gauge hierarchy prob- 
lem and, from a more phenomenological point of view, they would potentially be tested 
at the Large Hadron Collider (LHC). 

A particularly interesting scenario provided by theories in extra dimensions is the 
possibility of identifying the Higgs field with the internal component of a gauge field 
0. Several models of this sort have recently been considered, with [3] or without [U [3] 
supersymmetry, in various set-ups. Models invoking this scenario, given the origin of 
the Higgs field, are sometimes called models with gauge-Higgs unification. Electroweak 
symmetry breaking occurs radiatively in these theories, and is equivalent to a Wilson line 
symmetry breaking [Tj |H] , since the vacuum expectation value (VEV) of the Higgs field 
is proportional to the value of the Wilson line phase. In minimal 5D non-supersymmetric 
single-Higgs models, the Higgs potential is entirely radiatively generated and non-local in 
the extra dimensions. The models are highly predictive and the gauge hierarchy problem 
is solved. Matter can also be introduced in this set-up, with the possibility of naturally 
having a hierarchy of fermion masses jSJ. Although some problems have still to be solved 
before getting a completely realistic model, 1 theories with gauge-Higgs unification are 
definitely an interesting framework for physics beyond the SM. 

Aim of the present paper is to study in some detail the dynamics of 5D Wilson line 
phases at finite temperature, focusing in particular to the study of the Higgs potential in 
minimal S 1 /Z2 orbifold models with gauge-Higgs unification. By studying the one- loop 
effective potential, we find that the gauge symmetry broken by the Wilson lines is always 
restored at temperatures of order T ~ 1/L, where L is the length of S , in agreement 
with previous results jlUj . The behaviour of the effective potential at finite temperature is 
different from that usually expected in 4D or higher dimensional theories. This is best seen 
by looking at the dependence on the temperature of the mass of the field representing the 
order parameter of the phase transition. At high temperatures, it is known that the leading 
boson and fermion contributions to that mass are proportional in 4D to T 2 whereas 
in extra dimensional theories with n compact toroidal dimensions they are proportional 
to (LT) n T 2 ^2j, where L is the length of the extra circles, taken all equal, for simplicity. 
In our model the order parameter is given by the VEV of the Wilson line phase. Its non- 
local nature is responsible for a radically different behaviour at high temperatures: the 
fermion contribution is exponentially suppressed with the temperature, whereas the boson 
contribution has a leading term proportional to T/L, plus terms exponentially suppressed 
1 See however ref. [13 for a potentially realistic recent proposal. 
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with T. The mass of the Wilson line phase at zero temperature is radiatively induced and 
proportional to 1/L 2 ; when the latter is negative, a symmetry breaking is generated, but 
then at temperatures of order T ~ 1/L a phase transition occurs and the symmetry is 
restored. This will always happen, due to the presence of the bosonic contribution of the 
5D gauge fields. 2 The phase transition is analytically studied in a general framework, in 
presence of bulk bosons and fermions with generic couplings to the Wilson line phase. The 
transition is predicted to typically be of first order, with a strength which decreases as the 
fermion couplings increase. The appeareance of a first-order phase transition is essentially 
due to the presence of a term which is cubic in the Wilson line phase. The latter arises 
only from the contributions of massless 5D bosons. The phase transition takes place at 
temperatures of the same order of the compactification scale and thus well below the UV 
cut-off scale A of the model. 

Having established the generic presence of a first-order phase transition around T ~ 
1/L for Wilson line based models, we study some properties of the electroweak transition 
occurring in 5D models of gauge-Higgs unification. For definiteness, we focus to the 
particular model (and its generalizations) proposed in ref. 

Before summarizing our results, it might be useful to review very briefly the nature of 
the electroweak phase transition — and how it can reliably be studied — in the SM |14j . 
The latter properties are known to essentially depend on the value of the Higgs quartic 
coupling A, or equivalently on the Higgs mass Mh- For Mh < Myy, the transition is 
always of the first order, with a strength that is inversely proportional to the Higgs mass 
|15U16| . In this regime, the SM phase transition can reliably be computed in perturbation 
theory, although the resummation of an infinite class of certain IR divergent higher loop 
diagrams (so called "daisy" or "ring" diagrams) has to be performed. Also in the SM case, 
the presence of a first-order phase transition can be traced back to the presence of a term 
in the potential which is cubic in the Higgs field. Since the effective potential is generally 
a gauge-dependent quantity, care has to be taken in the gauge-fixing procedure JZj or 
a gauge-invariant formalism has to be advocated ^B]- Around the critical temperature, 
perturbation theory is less and less reliable as Mh approaches Mw- For Mh > Mw, 
perturbation theory breaks down and one has to consider non-perturbative methods, such 
as lattice computations. The latter seem to indicate that in this regime of Higgs masses, 
the SM does not actually have a phase transition but rather a crossover, since the Higgs 
and the confining phase are continuously connected |19j . 

One of the main motivations to study the detailed nature of the electroweak phase 

2 The restoration at high temperatures of a symmetry broken at zero temperature is not a universal 
feature. It has recently been observed, for instance, that in so called little Higgs theories the electroweak 
gauge symmetry remains broken even at high temperatures, at least up to the maximum temperature 
which can be studied in the effective theory I13| . In non-abelian theories, it can also happen that the 
Higgs and the confining phases are connected without the appearance of any phase transition. This is 
what lattice simulations predict for the SM, as we briefly remind below. 
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transition in the SM has been the realization that, under certain assumptions (namely 
a sufficiently strong CP violation and first order phase transition), it could satisfy the 
necessary requirements for baryogenesis [213] ■ Given the actual bounds on the value of the 
Higgs mass and amount of CP violation in the SM, it has definitively established that the 
SM phase transition is unable to give a successful baryogenesis. 3 

Let us now come back to analyze our results. In the model discussed in ref. 0, the 
Higgs mass is predicted and its value is typically Mu < My/ jl. In this regime, it is 
reasonable to compare our results with those obtained in the SM at one-loop level, for the 
same value of Higgs and top mass. 4 We find that the electroweak phase transition is of 
first order, with a strength that is inversely proportional to the Higgs mass, similarly to 
what happens in the SM (see figures . 

When 5D bulk fermions in very large representations of the gauge group are inserted, 
so that realistic values of the Higgs mass can be obtained, the first order phase transition 
becomes weaker and weaker (see figures On the other hand, one can get realistic 

values of the Higgs (and top) mass also by introducing large localized gauge kinetic terms, 
as discussed in ref. In this case the strength of the phase transition is considerably 
increased (see figures 113116(1 and could be strong enough to produce a reasonable baryo- 
genesis. A definite result can however be established only by a more careful study of the 
phase transition, which includes the study of sphalerons and of the possible sources of CP 
violation in our 5D set-up. 

As we have already mentioned, all our results, for any value of Mjj < My/ or Mu > 
My/ > are based on perturbative studies, obtained by computing the one-loop Higgs effective 
potential. In light of the importance of the correct gauge-fixing procedure, of the higher- 
loop daisy contributions and of the breakdown of perturbation theory for Mb > My/ in 
the SM, it is natural to ask if and to what extent one can trust our perturbative one-loop 
results. To this aim, we show that our 5D one-loop Higgs effective potential is considerably 
stable under radiative corrections. First of all, we establish that the 5D one-loop Higgs 
effective potential is gauge independent, by studying a class of background field gauges. 
This is expected, since the effective action typically fails to be gauge invariant away from 
its minima, but at tree-level in our model any VEV of the Higgs field is a good minimum, 
being absent any potential at this order. 

We give support to the idea that the phase transition can safely be studied in pertur- 
bation theory in Wilson line based models by computing the leading IR divergent daisy 
contributions to the Wilson line effective potential in a simple 5D model compactified on 
S . By a comparison between the naive one-loop thermal mass with the one-loop im- 

3 Proposals to improve the SM situation introducing effective dimension six (H^ H) 3 /M 2 operators 1211 
or strongiy coupied fermions 1221 have recentiy appeared. 

4 The top quark mass in the modei of ref. |5] is typically lower than the actual experimental one, 
M top ~ 174 GeV. 
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proved one — obtained by resumming an infinite class of daisy diagrams — we are able 
to establish that the leading daisy diagrams give a small correction up to temperatures 
of order T ~ 1/L (see figure 19 (a)) and start to give a relevant contribution only at 
temperatures of order T ~ 10/L or higher. We take these results as a strong evidence 
that perturbation theory in 5D models with gauge-Higgs unification is valid also around 
the critical temperature. The same comparison of one-loop naive and improved thermal 
masses is performed in the 4D theory obtained by a trivial dimensional reduction of the 
5D model, where only the Kaluza-Klein (KK) zero modes of the fields are retained. This 
truncation spoils the higher dimensional symmetries responsible for the finiteness of the 
Wilson line effective potential and indeed, as expected, daisy diagram contributions are 
larger (see figure 19 (b)). 

The paper is organized as follows. In section 2 we show the main features of one-loop 
effective potentials of Wilson line phases, namely their UV finiteness, gauge independence, 
high temperature behaviour and establish the generic presence of a first-order phase tran- 
sition. In section 3 we study the electroweak phase transition in the model of ref. [5| , which 
we briefly review in subsection 3.1. In section 4, the leading higher loop daisy diagrams in 
an U(l) model coupled with scalars and compactified on S 1 (and its corresponding model 
trivially reduced in 4D) are computed. Our conclusions are reported in section 5. Finally, 
we collect in an appendix all the formulae regarding the contribution to the Wilson line 
effective potential of bulk fields, with or without possible localized couplings. 

2 The Wilson line dynamics at finite temperature 

The form of the one-loop effective potential at T = for Wilson line phases on S 1 or 
S l /7j2 orbifold models is well-known [7JI2II]. The non-local nature of Wilson lines imply 
that their potential is completely UV finite, modulo irrelevant constant vacuum energy 
terms. This is the most important property of these potentials, which makes the idea of 
identifying the Higgs field with a Wilson line phase such an interesting proposal to solve 
the gauge hierarchy problem. We consider below only the effective potential for a single 
Wilson line phase. Effective potentials for multiple Wilson line phases, however, do not 
present any problem and are easily obtained by a generalization of that of a single Wilson 
line phase (see e.g. refs. [71 12*11 [2*5] ). 

In this class of models, gauge invariance severely constrains the couplings of the gauge 
field component A§ with other fields. In particular, any bulk field couples to A§ only 
through the minimal gauge coupling. In a suitable basis for the fields, each mode (i) couples 
diagonally with A$, giving rise to a simple shift of their KK mass (see e.g. appendix A of 
ref. [H] for details; all our conventions and notations are as in ref. [JJ]): = (n+q^a)/R, 
where n is the KK level of the state, q^> is the charge eigenvalue of the mode (i) and a 
(0 < a < 1) is the VEV of the Wilson line phase. It is related to the zero mode H of A§ 
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by the relation H = 2ia/(g4R), 5 where g\ is the 4D gauge coupling constant, related to 
the dimensionful 5D gauge coupling constant as follows: <?4 = g^/yL, with L = 2ttR. 

The effective potential at finite temperature, in the imaginary time formalism, is ob- 
tained from that at zero temperature by simply compactifying the Euclidean time direction 
on a circle of radius 1/(2ttT), where T is the temperature. Fields satisfying Bose/Fermi 
statistic must be taken to satisfy periodic/anti-periodic boundary conditions along the 
Euclidean time direction. The bosonic and fermionic contributions to the one-loop Wilson 
line effective potential at finite temperature are then easily evaluated. For each bosonic 
and fermionic degree of freedom with charge q and bulk mass term M, one has 



V(T,qa) 



(_)2r, T 



+00 



k,m= 



d 3 p 
(2tt)3 



log <^ p 2 + 



2ix{m + i])T 



+ M 2 (qa) 



(2.1) 



where M|(ga) = M 2 + (k + qa) 2 /R 2 are the masses of the states at KK level k, including 
a possible bulk mass term M, m are the Matsubara frequencies and rj = or 1/2 for 
bosons and fermions, respectively. It is important to recall that ea. (|2.1j) is valid not only 
for compactifications on S 1 but also on orbifolds, such as S 1 /^- In the latter case, 
the effect of the parity projection in the diagonal basis for the fields amounts simply to 
reducing the number of physical degrees of freedom with a given charge q with respect 
to the Wilson line phase. We collect in the appendix several equivalent ways in which 
one can compute the effective potential (|2.1|) . Notice that 60. 1)2.1(1 is valid also for gauge 
fields Am (M = 0, 1, 2, 3, 5); in this case one has to take into account a multiplicity factor 
equal to 3, the number of transverse polarization of a 5D gauge field. Moreover, the gauge 
contribution is gauge independent. We will discuss in more detail this important aspect 
in next subsection. 



2.1 Gauge independence of the effective potential 

The effective potential V(Hq) in standard spontaneously broken gauge theories, such as 
the SM, where Hq is the VEV of the Higgs field, is generally gauge-dependent This is 
easily seen at one-loop level, where the effective potential is obtained by evaluating a trace. 
In this case one can see, working for instance in an i?£-gauge, that the contributions of the 
would-be Goldstone bosons, ghosts and longitudinal gauge bosons, outside the minimum 
for the Higgs field, do not cancel anymore and leave a non-trivial ^-dependent, and thus 
gauge dependent, term. 

In gauge-Higgs unification models, where the Higgs field is identified with the KK zero 
mode of the gauge field A5, such problem does not occur, at least at one- loop level. This is 
expected since at tree level A§ is a modulus, namely any constant value of this field satisfies 

J The precise coefficient relating a and H is model dependent. For definiteness, we have taken here and 
in the following the one appearing in ref.|5]- 
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the equations of motion. In order to show more explicitly the gauge independence of the 
one-loop effective potential, we define in the following a class of gauge-fixing Lagrangians 
that are a sort of generalization of the 4D t'Hooft background gauges: 

1 



TV (D^ + iD b A 



5\2 



(2.2) 



where fj, = 0, 1, 2, 3 runs over the 4D directions and Dm = Qm + *<75[Am") •] is the covariant 
derivative in terms of the classical field configuration Am = Sm^oi/ {g^R)t, t being the 
direction in group space where the VEV is aligned. For £ = 1, the gauge-fixing term (|2.2() 
reduces to the usual 5D background field gauge commonly used in the literature, to derive 
the one-loop effective potential of Wilson line phases [7j; for £ = 0, we instead get the 4D 
Landau gauge d^A 11 = 0, whereas for £ — > oo, ea. Q2.2|) gives D$A 5 = 0; the latter condition 
is precisely the unitary gauge in which all would-be Goldstone bosons associated to the 
non-linearly realized gauge symmetries are decoupled, as we will see more explicitly in the 
following. From a 4D point of view, ea, ()2.2p is precisely the gauge-fixing term defining 
an i2g-gauge for the infinite gauge symmetry groups associated to all KK levels of the 5D 
gauge fields A^(x,y). The full 5D Lagrangian for the gauge fields is then 

1. 



C 



-Tr (F MN F MN ) + £ s + C 



gh > 



where C g h is the ghost Lagrangian associated to the gauge-fixing (|2.2|) : 



C 



-Tr[u>*{D fl D' 1 + 5 D° 



(2.3) 



(2.4) 



It is not difficult to derive the quadratic Lagrangian, in momentum space, for the gauge, 
scalar and ghost fields A^, A§ and u. As discussed before, the only non-trivial task is 
to derive the linear combination of fields that couples diagonally with A$, in which base 
the effect of the background is simply to shift the KK level of the mode. For definiteness 
(and because this will be the case we analyze mostly), we consider an SU(3) gauge theory 
on the orbifold 5 1 /Z 2 , where SU(3) is broken by the orbifold parity to SU(2)l x U(1)y 
and thus to U{\)em by the Wilson line phase a. The quadratic Lagrangian of this model 
reads, in momentum space, 



C 



quad. 



E E 

!=l,2n=- 

EE 

i=3 4 n=0 



A fi,n L ^tiv,n A v,n ' ^5,n ^55,71^5,71 ^ ^55,71^71 



+ 



A (i) A(0) A(i) , 4« a(0) A (i) +u (.)«a(0) (i) 
■ fi /i,n L1 fiu,n- fi i/,n ^ ^5,71^55,71^5,71 ^ w n L -^55,n UJ n 



(2.5) 



In eq. (|2.5j) . the fields (A^' , A$ , a/ 1 ' 2 )) represent the linear combination of fields that 



(3,4) 4 (3,4) w(3)4) 



couples diagonally to the Higgs field, with charge q = 1, 2. The fields (A^ ' , A5 
instead, do not couple to A5. The explicit form of these linear combinations of fields, in 
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terms of the usual SU(3) Gell-Mann field decompositions, can be found in Appendix A. 3 
of ref. [5]. Notice that due to the different Z2 parities of and ^5, we have A§ q = 



.4 



(i) 
^,0 



(i) 



0. The inverse propagators are 
(n + qa) 2 



A (3) 



p 2 + 



A 



(9) 

55, n 



P 2 + C 



# 2 

(n + qa) 2 
E 2 



+ 



-PfiPu 



9 = 0,1,2, 
9 = 0,1,2, 



(2.6) 



where p 2 = p^p^- It is clear from ea. ()2.6|) that when £ — > 00, all scalar fields A§ jn decouple 
from the theory, with the only exception of A^, namely the Higgs field. It is now easy to 
establish the independence on the parameter £ of the one-loop effective potential. Indeed, 
for each KK mode level n and charge q, the gauge, scalar and ghost contributions give a 
total factor 



A 



(?) 

55, n 



A (<?) detA (?) 



1/2 



VI[p 2 + 



(n + qaf 



(2.7) 



where the determinant in ea. (|2.7|) is meant to be taken only on the 4x4 polarization 
matrix of the gauge fields. Aside from the irrelevant factor, which is reabsorbed in the 
functional measure, all the non-trivial ^-dependence in eqs. ([2.6[) cancels out in eq. (|2.7[) . 
leaving an effective contribution of three scalar degrees of freedom with twist q. As stated 
before, this equals the total contribution of a 5D gauge field Am with the same twist. 

The gauge independence of the potential is trivially extended at finite temperature. 
Indeed, in the imaginary time formalism, ea. (|2.7|) still holds, the only modification being 
given by the Matsubara frequency modes, that are all integer-valued for gauge, scalar and 
ghost fields. 



2.2 High temperature behaviour 

The one-loop effective potential for Wilson line phases at finite temperature has been 
computed in ref. JHj (see also ref. |2S]), where it has been shown that at temperatures 
of order T ~ 1/L the symmetry is restored. Its high temperature behaviour is peculiar, 
as can be seen by studying the bosonic and fermionic contribution to the Higgs thermal 
mass at a = 0: 



M 2 H (T,a = 0) 



g A RY d 2 V 



(2.1 



a=0 



2 ) da 2 

For T S> the mass ()2.8(l can easily be extracted from ea. (|A,3|) . using the known 

asymptotic behaviour of Bessel functions for large values of their argument. For massless 
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bosonic and fermionic bulk fields (M = in ea. (|A.3|0 with twist q, one finds respectively 6 

,2 



(2.9) 
(2.10) 

where n b and n f denote the effective number of bosonic and fermionic degrees of freedom 



MUT,0)= - I ^(2V2vr 2 gV)^(LT)ie-^ 



{LT)i e -^ LT 
l + O 



1 

It 



of the field. The exponential terms appearing in eqs. (|2.9j) . (|2.10j) and the linear temper- 
ature dependence of the boson contribution are a characteristic feature of the non-local 
nature of the Higgs potential. They are different from the typical contributions expected 
for models in extra dimensions where one gets a thermal mass of order LT 3 |12| . or for 
the 4D Standard Model, in which the Higgs thermal mass, at high temperatures, goes like 
T 2 . 

The leading high temperature dependence of the one-loop thermal mass can be under- 
stood by studying the UV behaviour of the one-loop mass correction at zero temperature. 
Since this is in general quadratically divergent, power counting implies that its leading fi- 
nite temperature contribution will be of order T 2 . This counting can also be used in models 
in extra dimensions. Since at a given T, the effective number of KK modes contributing 
to the thermal ensemble is roughly given by LT, one would generally get a leading high 
temperature contribution to the thermal mass of order (LT)T 2 = LT 3 , in agreement with 
ref. |12) . For Wilson line potentials, however, the analysis is different, since they are finite 
at zero temperature, with a one-loop mass saturated by the compactification scale l/L 2 , 
times g 2 /(l6ir 2 ), a 4D loop factor. The non-local nature of the potential is such that all 
Matsubara mode contributions, but the zero mode one, are exponentially suppressed at 
high temperatures. The only relevant contribution arises then from bosons, which admit 
a zero Matsubara mode. Consequently, the thermal mass has a linear dependence on T, 
reproducing the leading term in ea. (|2.9|) . Since the effective number of KK modes con- 
tributing in the thermal ensemble is still LT, ea. (|2.9|) can alternatively be understood 
as the sum of the one-loop, zero temperature, mass corrections of the (LT) KK modes: 
(LT)/L 2 = T/L. The naive power counting argument, applied instead to a truncated sum 
of KK modes, would predict the wrong LT 3 dependence. In fact, a naive truncation of 
the KK sum would spoil the shift symmetry a — > a + 1 of the Wilson line phase, resulting 
also in a fake divergent zero temperature effective potential. 

The high T contribution to the Higgs thermal mass of massive bulk fields, with mass 

M, is further exponentially suppressed. In particular, for M 3> l/L, the linear dependence 

on T appearing in the first term in ea. (|2.9|) is suppressed by a factor e~ ML . The same is 

true for the whole effective potential, as is clear from ea. (|A.3|) . The effective potential is 

6 Recall that theories in extra dimensions are non-renormalizable and thus are effective theories valid 
up to an UV cut-off scale A 3> 1/R. For this reason, one cannot consider temperatures larger than (or of 
the same order of) A, because at these scales the theory becomes strongly coupled. 
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thus exponentially insensitive to UV physics and completely determined by the light 5D 
degrees of freedom. 

Given the high T behaviour of the bosonic and fermionic contributions to the Higgs 
mass in eqs. 1)2.9(1 and 1)2.10(1 . we can establish on quite general grounds that if a gauge 
symmetry breaking occurs at T = 0, namely M 2 (T = 0, a = 0) < 0, it will be restored 
at a sufficiently high temperature when M 2 (T,a = 0) > 0. Since Wilson line effective 
potentials are gauge invariant, V(— a) = V(a), and thus a = is always an extremum 
of the potential. Consequently, at sufficiently high T, a = turns to a minimum of the 
potential, at least locally. 

2.3 Properties of the phase transition 

An analytical study of the nature of the phase transition occurring in Wilson line based 
models is not totally straightforward, due to the non-local nature of the order parameter. 
The essential features of the transition, however, can quantitatively be studied by properly 
approximating the exact one-loop formulae for the potential appearing in the appendix. 

As we will show, the phase transition is typically of first order and is related to the 
presence of a cubic term in the bosonic contribution of the effective potential. This is very 
similar to what happens, for instance, in the SM, where a cubic term in the Higgs field 
appearing in the effective potential at finite temperature is responsible for the generation 
of a first-order phase transition, for sufficiently low values of the Higgs mass. The cubic 
term in the Wilson line phase a appears from a high temperature expansion of the bosonic 
contribution in, say, ea. ()A,3J) . For massless fields with M = 0, charge qs and for LT > 1, 
we can safely drop all Matsubara modes but the zero mode in the sum over m appearing 
in ea. ()A.3|) . For each degree of freedom, the resulting potential is 

V B (T, a) ~ ^ p . (2.11) 

k=l 

Interestingly enough, the potential 1)2.11(1 admits a very simple polynomial expression in 
a, which can be found by carefully Taylor expanding it in powers of a. Since a naive 
expansion in a would lead to the appearance of an infinite number of ill-defined sums, a 
regularization has to be taken to give a meaning to the resulting expression. We will use 
in the following the (^-function regularization, according to which £(— 2n) = J^fcLi = 
for all positive integers n. As a consequence, in the expansion of the cosine function in 
ea. ()2.11() . all terms but the first three, constant, quadratic and quartic in a, vanish. All 
the odd derivatives in the expansion vanish since £(— 2n + 1) is finite for n > 0, but care 
has to be taken for the third derivative. The latter in fact does not vanish, due to the 
relation 

lim > — i — '- = -. 2.12 

n=l 
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This is at the origin of the a 3 term in the expansion of the potential Q2.11J1 . It only 
arises for massless 5D bosonic fields, since for massive fields the sum over n is absolutely 
convergent, resulting in a trivially vanishing third derivative of the potential. 
Putting all terms together and recalling that £(0) = —1/2, one gets 



yp(T ' a) = -3l3 



— - (q B a) 2 + 2(q B af - (q B a) 



(2.13) 



The expressions 1)2.11(1 and (|2.13|) are identical for < a < l/q B . Since Vp(T,qpa = 
0) = Vp(T,qsa = 1), by periodically extending Vp to all values of a, one gets a complete 
identity between the two expressions. Indeed, it is straightforward to check that the 
Fourier series of the polynomial (|2.13|) agrees with ea. (|2.11() . Notice that the a 3 term 
is a non-analytic term, since it corresponds to a (H^ H) 3 / 2 operator, when written in a 
manifestly gauge-invariant form. 

Let us now study how the phase transition schematically occurs in a model with one 
massless 5D gauge boson and fermion, with charge respectively q B and qF with respect to 
a. We take qF > qp, because in this case the analysis is further simplified, as we will see. 
There is no need to specify, in this simple set-up, which is the underlying theory which 
gives rise to these fields and to the Wilson line. As a consequence, the analysis is general 
and it can be applied for any gauge group and compact space (S 1 or S 1 /Z2). 

At zero temperature, the effective potential is given by the first row in eq. ()A.5|) with 
M = 0. By summing the fermionic and bosonic contributions, we get 

3 ^ [8 cos(2irkq F a) - 3cos(2-irkq B a)] 

a) = ^ p • (2 ' 14) 

k=i 

The potential is unstable at a = and develops a minimum at a ~ l/(2qp), where the 
"Higgs mass" is approximately given by 

<^ 1^2^24^(3), (2.15) 

where, as rough approximation, we have neglected the bosonic contribution in ea. (|2.15|) . 
At finite temperature, as we have already discussed, the fermion contribution is less and 
less relevant, resulting in a phase transition for T ~ 1/L. The transition is more efficiently 
studied by starting from temperatures slightly above the critical one, when TL > 1. 
In this regime, it is a good approximation to use the high temperature expansion for 
the potential. The bosonic contribution has been already computed and gives rise to 
the polynomial potential (|2.13J) multiplied by a multiplicity factor of 3. The fermionic 
contribution, unfortunately, does not admit a simple finite polynomial expansion which 
accurately approximates the exact one-loop result, like the bosonic case. As such, a rough 
approximation will be used to put it in a simple polynomial form. By expanding the 
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modified Bessel functions and retaining only the modes m = 0,-1 and k = 1 in ea. l[A.3[) . 
we roughly get 

Vf(T, a) ~ ^(LTf' 2 er vLT cos(2vrg F a) . (2.16) 

No a? term arises from the fermion contribution. If we expand the cosine in ea. (|2.16|) and 
retain only terms up to a 4 we see that the effective potential for LT > 1, resulting by 
summing Vb and Vf, has the schematic form 7 

L A 

-~V(T, qa) ~ a(x) a 2 - b(x) a 3 + c(x) a A , (2.17) 

7T^ 



where x = LT and 



o(x) = - S^v^x 5 / 2 --™ 



b(x) = 2xq% , (2.18) 
c(x) =4x + ^4v / 27r 2 x 5 / 2 e-™. 

Eg. (|2. 17(1 is valid for < a < l/(2qp), which is the relevant range in a for the study of the 
phase transition. The analysis of the latter for potentials like ea. (|2,17|) is standard (see e.g. 
ref. |27|). For x >> 1, the potential admits only a minimum at a = 0. As the temperature 
decreases, an inflection point appears at T = T\, below which a non-trivial minimum 
and maximum appear. The critical temperature Tq is defined as the temperature when 
V(Tc, qot m in{Tc)) = 0, namely when the non-trivial minimum is as deep as the minimum 
at zero. This is given by the largest root of the equation b 2 (xc) — ^a(xc)c{xc) = which, 
for qp > qs, is nearly equivalent to the vanishing of the term a(x) in ea. Q2.17[) . The critical 
temperature has only a logarithmic dependence on the charges qs and qF and is of order 
1/L, as mentioned. At T = Tc, we get 

a (xr) ~ b{Xc) ~ (2 19) 

a mm (x c ) - 2c{xc) - ^ . (2.19J 

Below the critical temperature Tq, a m i n (T) becomes the new global minimum of the 
potential. The value \H(Tc)\/Tc = 2a m i n (Tc) / (g^RTc) is one of the relevant parameters 
to study, if one wants to get baryogenesis at the electroweak phase transition (see e.g. 
ref. EE])- One necessary (but not sufficient) requirement is that \H(Tc)\/Tc > 1, otherwise 
sphalerons at T < Tc would wash out any previously generated baryon asymmetry [201 - 8 ^ 
is also a good parameter to measure the strength of the first-order phase transition. Since 

7 As we said, this is generally a rough approximation but it is enough for our analytical estimates. All 
our results in the following section are instead based on numerical analysis which take into account the 
exact form of the one-loop boson and fermion contribution to the potential, as reported in the appendix. 

8 The bound \H(Tc)\/Tc > 1 is derived by an analysis of sphaleron dynamics in the SM. We have not 
repeated this analysis for Wilson lines in 5D and thus deviations from the bound H(Tc)/Tc > 1 could be 
present. 
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the Higgs mass ()2.15j) is proportional to qp, whereas a m i n {xc) ~ qs/qp anc ^ ~ l/A 
we conclude that 

\H(T C )\ q B 
Tc q F 

The strength of the first-order phase transition is inversely proportional to q\ and hence 
to the value of the Higgs mass. 

In presence of several massless bosonic and fermionic fields, with charges qs,i and qp,i, 
the potential can still be put in the form ()2.17j) . provided that one substitutes in eas. (|2.18|) 
q k B — > nB,iq% j and q F — > Yli n F,iqpi, where riB(F),i are the multiplicities of the bosons 
and fermions with charges qs(F),i = 2,3,4). 

For more generic field configurations, such as massive bulk fermions, localized fields or 
localized gauge kinetic terms, the analytical study of the phase transition is more involved 
and one has to rely on numerical methods to safely establish its nature. However, since 
the coefficient b(x) multiplying the a 3 term in ea. (|2.17|) will still be non- vanishing (unless 
the model has no massless bosons coupled to the Wilson line phase), we expect that the 
phase transition will be of first order. We show in next section that this is what generally 
happens for the model of ref. jS], by numerically studying the one- loop Higgs potential in 
a more complicated setting. We will also see that the analytical studies performed in this 
subsection are still approximately valid. 



3 The electroweak phase transition in gauge-Higgs models 

In this section, we want to establish the nature of the electroweak phase transition in 
models with gauge-Higgs unification based on 5D orbifolds. We will consider in particular 
the model constructed in ref. [H] and study if and how the electroweak phase transition 
occurs in it and some of its extensions. 



3.1 Review of the model 

The model of ref. is based on a 5D gauge theory with gauge group G = SU (3) c x SU (3) w 
on an S 1 j'Li orbifold. 9 The Z2 orbifold projection is embedded non-trivially only in the 
electroweak SU(3) W gauge group, by means of the matrix 

/— 1 \ 

p = e mXa = 0-1 , (3.1) 
V 1/ 

where A a are the SU(3) Gell-Mann matrices, normalized as TrA a A& = 25 a b- The twist 

(|3.1|) breaks the electroweak gauge group in 4D to SU(2) x U(l). The massless 4D fields 

9 As explained in detail in ref. [5], an extra 17(1)' gauge field has to be introduced in 5D, in order to get 
the correct weak mixing angle. This additional gauge field does not couple with the Higgs field, and is thus 
irrelevant in all the considerations that will follow. For this reason it will be neglected in the following. 
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are the gauge bosons in the adjoint of SU{2) and a charged scalar doublet H, the 
Higgs field, arising from the internal components A\ of the gauge field. A VEV for A\ 
induces an additional spontaneous symmetry breaking to U(1)em- Following ref. j^j, we 
can take 



corresponding to a purely imaginary Higgs VEV (H) = 2ia/{g^R), where a is a Wilson 



The exact symmetry a — > a + 1 prevents the appearance of any local potential for the 
Higgs field, at any order in perturbation theory. The only allowed potential terms for the 
Higgs are non-local in the extra dimension (Wilson lines) and thus free of any divergence 
|29j . Despite the model is a non-abelian gauge theory in 5D, and thus non-renormalizable, 
the Higgs potential is insensitive to UV physics and thus calculable. The gauge hierarchy 
problem is automatically solved in these models. 

Introducing matter fields in this set-up is a non-trivial task. One possibility is to 
introduce massive 5D bulk fermions and massless localized chiral fermions, with a mixing 
between them, so that the matter fields are identified as the lowest KK mass eigenstates. 
In this way, a realistic pattern of Yukawa couplings can be obtained [3]- 10 

The electroweak phase transition occurs in these models because the radiatively in- 
duced one-loop Higgs potential, at zero temperature, is unstable at a = 0, in presence 
of matter fields. In the broken phase, the mass of the matter fields is exponentially sup- 
pressed in terms of the bulk fermion masses Mj. Since the contribution to the Higgs 
potential of massive 5D bulk fermions is also exponentially suppressed in Mj, effectively 
the electroweak phase transition is driven by the lightest bulk fermion, the one giving rise, 
through its mixing with the corresponding localized chiral fermion, to the top quark and 
top Yukawa coupling. In this sense, the electroweak phase transition is driven by the top 
quark. We hence focus in the following to the study of the contribution to the Higgs effec- 
tive potential at finite temperature, given by the gauge fields and by a single bulk fermion, 
coupled to the boundary fermions which give rise to the top quark. Subsequently, in order 
to get more realistic values for the Higgs and top mass, as well as for the compactification 
scale 1/R, we will consider some extensions of this minimal model, as already suggested 
in ref. [Sj. 

The bulk 5D fermion that has the correct quantum number to couple with the (conju- 
gate) top quark is the symmetric representation 6 of SU (3),,,. 11 More precisely, the matter 

10 It should be said, however, that the construction of an UV completion leading to this fermion spectrum 
might be a non-trivial task. Indeed, the latter does not satisfy the requirements recently found in ref. |3Uf 
that allow to consider this model as a degeneration limit of a well-defined 6D model on T 2 /Z2. We thank 
A. Wulzer for this observation. 

"Thanks to the extra 17(1)' symmetry, it is actually possible to couple the top quark with larger 
representations of SU(3) W , resulting in this way to a larger value for the top mass. We will not pursue 




(3.2) 



line phase, < a < 1, and 54 = gs/VX is the 4D weak coupling constant. 
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fermion of this basic construction is as follows. We introduce a couple of bulk fermions VP 
and \F with opposite Z2 parities, in the representation (3,6) of SU(3) C x SU(3) W . At the 
orbifold fixed points, we have a left-handed doublet Ql = (tL,bi) T and a right-handed 
fermion singlet tR of SU(2) x U(l). They are located respectively at y = and y = L/2, 
the two boundaries of the segment S 1 /Z2. The parity assignments for the bulk fermions 
allow for a bulk mass term M mixing and ^f, as well as boundary couplings e\p with 
mass dimension 1/2 mixing the bulk fermion ^ to the boundary fermions Ql and i_R. The 
matter Lagrangian reads 



mat 



* 5 * + * + ( + h.c. 



+ <%) Q l 4 Ql+ (eiQ^V + h.c. 



+ % " *fl 04 t* + ( e 2 t£* + h.c. ) , (3.3) 



2 

where ip and x are the doublet and singlet SU(2) components of the bulk fermion fi± 
and 1)5 denote the 4Z? and 5.D covariant derivatives so that 1)5 = l^ + ^Ds. The exact 
spectrum of the bulk-boundary fermion system defined by the Lagrangian (|3.3j) can be 
determined by solving a complicated trascendental equation. The lighest state is however 
massless, as long as a = 0. For the mass M top of the lightest state, that we identify 

with the top quark, is dynamically generated by non-local Wilson-lines operators, whose 
coefficients are exponentially suppressed by the bulk mass term M. The W mass equals 
My/ = a/R, whereas the Higgs mass is radiatively induced and equals 



M 2 H {a r 



2 «2 



8 Z V 



( 94R 
V 2 J da 2 



(3.4) 



with V(a) the zero temperature Higgs effective potential and a m i n its minimum with the 
lowest value of a. This is typically the global minimum of the potential in the broken 
phase and is also the first one appearing during the electroweak phase transition. Even 
when a m i n is only a local minimum, semiclassical analysis have shown that its tunneling 
decay rate is negligibly small, leading to a practically stable minimum |24j . 

Finally, let us comment on the range of validity of this theory as an effective field theory. 
Using nai've Dimensional Analysis (NDA), we can estimate the UV cut-off scale A to be 
given by the scale in which the 5D gauge couplings become strong. Using respectively the 
strong SU(3) C and weak SU (3) w coupling constants, one gets A c ~ W/R and A^ ~ 100/i? 
0. 



3.2 Results 

We report below our study of the electroweak phase transition in the model of ref. |S] 
and its generalizations. The results, including all the figures, have been obtained by a 
this possibility. 
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numerical computation of the exact one-loop relations reported in the appendix. 
3.2.1 Only bulk fields 

This is the simplest set-up, in which we consider gauge fields and a couple of bulk fields, 
with opposite Z 2 parities, in the representation (6,3) of SU(3) W x SU(3) C . No boundary 
fermions are introduced. 

500 ■ : 




250 f, 

10 15 20 25 30 

M„ (GeV) 



Figure 1: Compactification radius as function of the Higgs mass (only bulk fields). 



In this theory, a first order phase transition appears for A < 2.1, where A = ML/2 is 
the properly rescaled bulk fermion mass, while outside this range no symmetry breaking 
occurs. The critical temperature of the transition is Tq ~ 1/L. We report in figures ^ 
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(a) 




20 25 
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(b) 



Figure 2: (a) Critical temperature and (b) position of the minimum at T = and at 
T = Tc, as function of the Higgs mass (only bulk fields). 



and |21 how the compactification radius, the critical temperature and the minimum of the 
potential (at T = and at T = Tc) depend on the Higgs mass, whose value depends on A. 
In figure |21 (a) , we plot \H{Tc)\/Tq as function of M#. The phase transition is strongly 
first order, as expected for such low values of Mjj- In figure 0(b), we plot the latent heat 
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as function of Mjj (normalized in units of the critical temperature), whose expression is 
given by 

(3.5, 



AQ=T—V(T,a min ) 







T=T C 



T—V(T,a = 0) 



T=T C 



For simplicity, we normalize the effective potential so that V(a = 0, T) = 0, in which case 
the second term in ea. (j3.5j) vanishes. Finally, in figure 0] we show the shape of the effective 
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Figure 4: (a) Effective potential (for Mh = 25 GeV) at various temperatures, (b) Detail 
near the phase transition (only bulk fields). 



potential, for Mjj = 25 GeV, at various temperatures. 

For high values of Mh (which correspond to low values of the bulk fermion mass A), 
the behaviour of a m i n (Tc) and \H(Tc)\/Tc in figures [21 (b) and 0(a) can be understood 
by applying the considerations of subsection 2.3. Around the critical temperature the 
potential 1)2.17(1 has only a mild dependence on A, at least for A < 1.5, since as first 
approximation it implies only the shift ttLT — > \/ 4\ 2 + (irLT) 2 in ea. (|2.16l) . Hence, 
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although the Higgs mass at T = is sensibly dependent on A, a m i n and \H(Tq)\/Tc, as 
computed in eqs. (|2.19|) and (|2.20j) . are essentially independent of Mh for low values of A 
(high values of Mh)- 

For larger values of A, both Tc and Mh starts to exponentially depend on A, resulting 
in a linear dependence of Tq on Mh- On the contrary, a m in{Tc) remains still constant. 
These behaviours are roughly reproduced by the left part of figures 0(b) and El (a). 

3.2.2 Inclusion of boundary fermions 

We now add the boundary fermions, namely the doublet Ql and the singlet tn, as in 
ea. (|3.3|) . The relevant formulae for the effective potential in presence of boundary fields 
are found in the appendix, eas. (|A.7|) . Since the potential does not significantly depend on 
the bulk-to-boundary coupling e^, the crucial dependence being given by the bulk fermion 
mass A, we have set £2 = ei and fixed the overall scale by the requirement of having 
Mtop = 45GeV. The effective potential is studied for 0.85 < A < 1.85. Again, the 
electroweak phase transition is of first order, with a critical temperature of order 1/L (see 
Figs. 15m 




1 2 50[ , , , , ^ 
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(a) (b) 

Figure 5: (a) Values of the e\ coupling and (b) compactification radius as function of the 
Higgs mass (inclusion of boundary fermions). 

For comparison, in figures|Hl(a) and0 we plotted the critical temperature, \H(Tc)\/Tc 
and the latent heat as function of the Higgs mass for both the 5D model and the SM, 
with M top = 45 GeV. For the SM case, we plot both the naive one-loop {SMjj nres ) and 
one-loop improved {SM£> a i sy ) results, the latter obtained by resumming the leading daisy 
diagrams ^3E3- 12 For such low values of the Higgs mass, 10 -j- 20 GeV, the SM effective 

12 We have fixed the Landau gauge, which turns out to be a good gauge choice. Indeed, results obtained 
in this gauge are very similar to the ones obtained in gauge-invariant approaches to the electroweak phase 
transition, such as the gauge-invariant effective potential 1181 and lattice computations |31| . We also used 
a cut-off renormalization scheme choosing the counterterms so that the position of the zero temperature 
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Figure 6: (a) Critical temperature and (b) position of the minimum at T = and at 
T = Tq, as function of the Higgs mass (inclusion of boundary fermions). 




Figure 7: (a) Phase transition strength and (b) latent heat as function of the Higgs mass 
(inclusion of boundary fermions) . 



potential is well described by its one-loop approximation. 

Even when boundary fermions are inserted, the behaviour of a m i n (Tc) and \H(Tc)\/Tc 
reported in Figures (b) and [7| (a) is very similar to that found in the previous case 
with only bulk fermions. This is explained by recalling, as already mentioned, that the 
potential has a small dependence on the bulk-to-boundary couplings e,,. All the analytical 
considerations performed at the end of section 3.2.1 are then approximately valid also in 
this case. 

minimum of the effective potential and the Higgs mass do not change with respect to their tree level values 
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Figure 8: (a) Effective potential (for Mh = 18 GeV) at various temperatures, (b) Detail 
near the phase transition (inclusion of boundary fermions). 



3.2.3 Theory with high rank bulk fermions 

One of the main problems of models based on gauge-Higgs unification is the predicted 
too low value for the Higgs mass. A possibility to solve this problem is given by the 
introduction of additional bulk fermions, which do not couple to the localized matter fields, 
in high rank representations of SU (3) w ■ We consider in the following the electroweak 
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Figure 9: (a) Value of the e± coupling and (b) compactification radius, as function of the 
Higgs mass (high rank bulk fermions). 



phase transition that one obtains in a 5D model where, in addition to the fields introduced 
in section 3.2.2, one considers a massive bulk fermion in the completely symmetric rank 
8 representation of SU(3) W and singlet of SU(3) C . The mass of this fermion has been 
fixed to Xhr = 2.15. The bulk-to-boundary couplings q have been fixed requiring that 
Mtop = 45 GeV, with €2 = 0.5 ei. The effective potential is studied for 1.6 < A < 2.4. 
For A > 1.79 (which corresponds to Mjj > 125 GeV), a m i n is the global minimum of 
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the potential, which turns to a local minimum for smaller values of A. Our results are 
summarized in figures 191121 The Higgs mass can now reach realistic values (although 
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Figure 10: (a) Critical temperature and (b) position of the minimum at T = and at 
T = Tci as function of the Higgs mass (high rank bulk fermions). 
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Figure 11: (a) Phase transition strength and (b) latent heat as function of the Higgs mass 
(high rank bulk fermions). 



the top mass is still too low) . As is clear from figure ^2 (a) , a very weak first order phase 
transition is found in this model with Tq ~ 1/L. 

This behaviour can be understood by means of our general analysis of section 2.3. 
The fermionic potential is dominated by the high rank fermion whose mass is fixed. Since 
Mh varies only with A, the small dependence on Mjj in figures (b) and ^2 is easily 
understood. The decreasing of the strength of the phase transition with the rank is clear 
by noting that 

a mm (T c )~|^4, (3.6) 
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where i runs over all the fermion components with charge qi and multiplicity 7ij, arising 
from the decomposition of the rank r bulk fermion. 

Fermions in too high-rank representations can not be considered, since the value of 
the UV cut-off A w rapidly decreases as r increases, restricting the range of validity of the 
effective theory [H]. No comparison with the SM is given, since for such values of the Higgs 
mass, perturbation theory breaks down close to the critical temperature. 
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Figure 12: (a) Effective potential (for Mh = 135 GeV) at various temperatures, (b) Detail 
near the phase transition. The effective potential is multiplied by 10 with respect to the 
(a) diagram (high rank bulk fermions). 



3.2.4 Theory with localized gauge kinetic terms 

An interesting possibility that is allowed by symmetries in our 5D model is the introduction 
of localized gauge kinetic terms. These have a significant effect on the phenomenology of 
models based on gauge-Higgs unification. In particular, they represent another option, in 
alternative to the introduction of fermions in high rank representations, to get realistic 
values for the Higgs mass. Interestingly enough, their presence also allow to get acceptable 
values for the top mass. As a drawback, they break the custodial symmetry, leading to 
too large deviations for the p parameter. 

We introduce localized gauge kinetic terms at only one fixed point of <S' 1 /Z2, since 
this is the most phenomenological convenient option. In the notation of ref. [5], we take 
ci = c = 6 and C2 = 0. 13 The gauge contribution to the Higgs effective potential, in 
presence of localized kinetic terms, is reported in ea. (|A.ll|) of the appendix. The fermions 
are as in section 13.2.21 The bulk-to-boundary couplings £j are chosen so that the top 
quark mass is 110 GeV, with €2 = e\. The effective potential is studied for 0.8 < A < 1.9. 
Interestingly enough, despite the high value of the Higgs mass obtained in this case, 

13 When localized gauge kinetic terms are inserted, the cut-off A,„ — > A w /c, so the value chosen for c is 
well below the limit of validity of the effective theory. 
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Figure 13: (a) Value of the t\ coupling and (b) compactification radius as function of the 
Higgs mass (localized gauge kinetic terms). 



110 4- 170 GeV (see figures IT3*1 and HI)) ■ the phase transition is still of first order and is 
also moderately strong, with \H(Tc)\/Tc ~ 0.7 (see figure ITol (a)). Since this is of order 
one, it might be strong enough to avoid that sphalerons in the broken phase wash out the 
previously generated baryon asymmetry (see footnote |HJ) . The critical temperature in this 
model is ~ 2/L (see figures H3*l (b) and ITU (a)). 
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Figure 14: (a) Critical temperature and (b) position of the minimum at T = and at 
T = Tci as function of the Higgs mass (localized gauge kinetic terms). 



The behaviour of the phase transition for large values of c is analytically harder to be 
studied than the previous cases, since the localized gauge field contribution (|A.11|) plays 
a crucial role. By expanding the total gauge field contribution, we have been able to 
establish that the first two terms of the potential in an a expansion, namely the quadratic 
and cubic term in eq. Q2.17jl . scales like a(x) ~ 1/c and b(x) ~ c~ 3 / 2 . We have not been 
able to find the asymptotic behaviour for large c of the quartic coupling, which seems to 
be negative. Neglecting the quartic coupling arising from the bosonic potential, we find 
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that the critical temperature depends only logarithmically on c and awn(^b) ~ Va/c- 
Since 55 ~ g^\fLc, we get H(Tc)/Tc ~ 1/c. Our numerical analysis confirm the decrease 
of the phase transition strenght with c, but numerically we find \H(Tc)\/Tc ~ l/c 1 / 2 " 1 , 
which indicates that the bosonic quartic coupling cannot be totally neglected and result 
in a slightly stronger first-order phase transition. 
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Figure 15: (a) Phase transition strength and (b) latent heat as function of the Higgs mass 
(localized gauge kinetic terms). 
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Figure 16: (a) Effective potential (for Mh = 135 GeV) at various temperatures, (b) Detail 
near the phase transition (localized gauge kinetic terms). 



4 Estimate of higher-loop corrections 

All the analysis performed in this paper about the electroweak phase transition is based 
on the computation of a one-loop effective potential. Aim of this section is to estimate 
the leading higher loop contributions and determine the range of validity of our one-loop 
analysis. 
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It is known that the most important higher loop contributions to the effective potential 
in 4D theories generally arise from certain IR divergent diagrams, known as "daisy" or 
"ring" diagrams (see figure fTTj) . 





m+2 



(b) 



Figure 17: Examples of higher loop daisy diagrams contributing to the mass of the Wilson 
line phase. Solid and wavy lines represent respectively the charged scalar and the Wilson 
line phase. For simplicity, we have omitted the arrows in the charged scalar propagators. 



As shown in ref . jS2] ; the latter are all easily resummed by adding and subtracting to the 
Lagrangian the one-loop (thermal) mass correction for any field present in the Lagrangian. 
One then defines a new propagator by adding the thermal mass to the tree level one, and 
consider the other (opposite) term as a one-loop mass insertion. Such procedure has 
been performed both at one JH] and two-loop level |16j in the SM. Interestingly enough, 
the two-loop improved SM effective potential turns out to give results very similar to 
that computed by numerical results based on lattice simulations even for values of Mh 
close to (but less than) Myy, in which case perturbation theory is expected to be a bad 
approximation . 

In the following, we will give an estimate of the contribution of daisy diagrams for 
models where the order parameter is a Wilson line phase, as in theories with gauge-Higgs 
unification. For simplicity, we focus on a particularly simple model, scalar QED in 5 
dimensions, compactified on a circle S 1 of length L = 2irR. u The La grangian is 

C = {D M <S>) ] D M <S> - ^-(&<S>) 2 - \F MN F MN , (4.1) 

where $ is a complex scalar field with periodic boundary conditions: Q(2ttR) = <£(()). 
Our aim is to study the effective potential for the Wilson line phase a = A^g^R, in the 
approximation in which the dimensionless 4D gauge coupling 54 <C A, with A = A5/L. 15 

14 We take S 1 and not S 1 /Z2 as compact space because in the former case we can study an abelian gauge 
field model, whereas in the latter case the Z2 orbifold parity forbids the presence of abelian Wilson line 
phases and would thus oblige us to consider non-abelian gauge fields. 

15 By performing a non-single valued (on S 1 ) gauge transformation, we could alternatively reabsorb the 
Wilson line phase a into the twist of the scalar field: &(2nR) = e 2 ' l7T °'^(0). 
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At one-loop order, the effective potential V(a, T) is given by one of the relations reported 
in the appendix, eqs. (|A.3| ). (|A.5|) or (|A.6j) . with d = 3, M = 0, r/ = and a factor of 
2 taking into account the two scalar degrees of freedom. At higher order, the leading 
contributions arise from the daisy diagrams drawn in figure El Since g& <C A, we can 
consistently neglect the daisy diagrams involving one-loop mass corrections induced by 
the gauge field and consider only the one arising from the quartic coupling A (see figure 
18). 

Let us denote by /„ the (n + l)-loop daisy diagrams obtained by summing all possible n 
one-loop mass insertions at zero external 5D momentum and define an effective parameter 
7(n) = I n /I n -i; if there exists some n for which j(n) ~ 1 or higher, perturbation theory 
breaks down and all daisy diagrams must be resummed. It is not difficult to compute 
7(n) for n> 1. In this limit, one has to compute only the diagrams (b) in figure IT71 
since they have a combinatorial factor proportional to n, whereas the diagrams (a) have 
a constant combinatorial factor. Moreover, the n one-loop mass insertions, namely the 
"small" bubbles in red in figure IT71 decouple from the computation and one has only to 
compute the finite "big" bubble, the one in black in figure fT7l (b). We get 



lin) 



LMi^i oop (T,a) 
2vr 



E 



(k + a) 2 + (LT) 2 m 2 



E 



m,k 



+ a) 2 + (LT) 2 m 2 



-2n + l 



(4.2) 



In eq. ([4.2j) . M 2 _ loop (T,a) is the one-loop finite thermal mass correction for the field $ 
which, for A ^> g, is given only by the diagram in figure IT8l and reads 

(LT)T 2 X^ ^ cos(27rafc) 

m 2 + (LT) 2 k 2 



M 



l—loop 



(T,a) 



4tt 2 



3/2 ' 



(4.3) 



The mass correction (|4.3|) is defined in a regularization in which only the finite terms 




Figure 18: The one-loop mass correction M 2 _ loop to the charged scalar <£, induced by the 
A($t(j))2 COU pii n g. 



arising from the compactification are left, neglecting the UV divergence taken to be equal 
to that obtained in non-compact 5D space. This is technically achieved by dropping the 
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Poisson resummed k = term in ea. (|4.3j) . We immediately recognize from ea. (|4.2|) that 
for a ~ 0, the leading contribution to I n is given by the KK and Matsubara zero modes 
k = m = 0, which give rise to |7(n)| ~ I /a 2 . The very low values of the Wilson line phase 
a are the most affected by daisy diagrams, as is obvious by the fact that IR divergences 
in the diagrams I n appear precisely in the limit a — > 0. A more transparent estimate of 
the relevance of resumming daisy diagrams can be performed in the two extreme limits 
when LT <C 1 and LT S> 1. In these two limits, and for sufficiently high n, ^{n) ~ 7 is 
independent of n. For LT <C 1, we can approximate the double sum in ea. (|4.2|) by taking 
k = and substituting the sum over the Matsubara modes with an integral. Viceversa, 
in ea. Q4.3JI . we can take rh = 0, retain only the sum over the modes k and fix a = 0. 16 In 
this way, we get, for LT <C 1, 

7 - 4vW " M f ree ' 1 • ) 

where M 2 ree = (a/R) 2 is the tree- level mass of the lowest KK mode for the field <1>. Daisy 
diagrams start to be relevant when the thermal mass is of the same order (or higher) than 
the tree-level one. 

For very high temperatures, LT 3> 1, the double sums in ea. (|4.2|) are well approximated 
by the leading term k = m = 0, whereas in the mass correction 1)4.3(1 we can approximate 
the sum over rh with an integral (setting again a = 0) and sum over all k. One gets, for 
LT » 1, 

XLT 1 IA * 

^"48^- (45) 

As eq. ([4.4[) . eq. (|4.5|) can be written as the ratio of the thermal mass at high temperatures 

over the tree-level mass M 2 ree . In fact, as a rough estimate, obtained by taking k = m = 

in the two double sums in ea. (|4.2|) . one sees that at any temperature 7 ~ — M 2 _ t oop /M 2 ree . 

Once established the necessity of considering higher-loop daisy diagrams for sufficiently 

low values of the Wilson line phase a, let us resum all these diagrams by following the 

procedure outlined in ref. |S2]- Since 7 < 0, one can expect some cancellation between daisy 

diagrams of different order, resulting in a decrease of the effect of the daisy resummation 

procedure. This is indeed what happens, as we will see. As briefly discussed before, 

along the lines of ref. |S2]) the resulting one- loop improved effective potential Vi mp (T, a) is 

obtained by one of the relations reported in the appendix, eqs. l)A.3| ). (|A.5|) or (|A.6|) . but 

16 Notice that for sufficiently high values of a, a > 0.2, Ml_ loop (a,T) becomes negative and a non- 
vanishing VEV for the field $ is induced, which complicates the above analysis. For this reason, we 
analyze here only the region of low values of a, which are anyway the ones most affected by the daisy 
diagrams. 
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by taking M = Mi 

— ioop(T,0) instead of M — 0. Prom ea. ()A.6| ). for instance, one gets 



rr lT v 4LT 5 ^ ~ cos(27rfcq) / A/i^ oop (T, 0) /~ 2 ^ r ^ 2 \, 

(27r) / ^ ^ [^2 + (LTfc)2] 5 V T V y 

(4.6) 

A good estimate of the higher-order daisy contributions at any temperature is per- 
formed by comparing the Wilson line thermal mass obtained from Vi mp (T, a) with that 
derived from V(T,a). We denote them respectively M 2 mp (T, a) and M 2 (T,a) (not to be 
confused with M 2 _ loop (T,a) in ea. (|4.3|) . which is the thermal mass of the field $). We 
plot in figure 19 (a) Mf (T, a) and M 2 (T, a) as function of a for A ~ 1 and for LT ~ 1. 
As can be seen, the two behaviours are very similar, with deviations which are at most of 
about 10% for a = 0. As expected, the maximal deviation is found for values of a ~ 0. 
We find that at low temperatures (LT ~ 1/10), the maximal deviations (at a = 0) are of 
order of a few %, whereas at high temperatures (LT ~ 10) they are of order 30%. Notice 
that, although this model does not have a phase transition in the Wilson line parameter, 
LT ~ 1 is in the range of the critical temperature where the phase transition occurs in 
such kind of models, as we have explicitly seen in sections 2 and 3. From our analysis we 
thus find that the leading higher-loop daisy diagrams give small corrections to the effective 
potential. This can be taken as an evidence of the fact that perturbation theory in our 5D 
models is still valid at T ~ Tq. On the contrary, for very high temperatures (LT >> 10), 
the one-loop improved potential starts to significantly differ from the naive one-loop po- 
tential. This is particularly evident from the equivalent expression of the potential given 
in ea. (|A.3|) . At very high temperatures, the naive one-loop potential is entirely dominated 
by the Matsubara zero mode m = 0, proportional to T, whereas the one-loop improved 
potential, due to the non- vanishing argument in the function B2 appearing in ea. ()A.3|) . 
is exponentially suppressed in LM ~ \fLST. The discrepancy that one finds between the 
two potentials starts however to be very large for temperatures which are well above the 
cut-off scale A. As discussed in footnote U3 it is meaningless to consider such range of 
temperatures in an effective field theory approach. 

It is useful to compare the contribution of daisy diagrams for the 5D U(l) model above 
with that of the 4D model arising by trivial dimensional reduction, where all the massive 
KK modes are neglected. By denoting with a and 4> the properly normalized zero mode 
components of A5 and <!>, one gets the following 4D Lagrangian: 

C = (D^D»4> + i(fya)(0"a) - g 2 a 2 - j(tfcf>) 2 - . (4.7) 

The trivial dimensional reduction spoils the 5D Wilson line nature of a, breaking the shift 
symmetry a — > a, + 1. Let us study the effective potential V(a) for the VEV of the field 
a ((a) = a/(gR)). The contribution of <E> to V(a) is the standard one for a massive scalar 
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Figure 19: (a) Thermal masses in 5D, from the one-loop naive potential (red solid line) 
and the one-loop improved one (blue dashed line), (b) Thermal masses in the reduced 4D 
model, from the one-loop naive potential (red solid line) and the one-loop improved one 
(blue dashed line). In both cases the normalization is such that M 2 (a = 0, T) = 1. 

field at finite temperature One has, neglecting terms independent of a, 



In eq. (|4.8|) . V(0, a) is the zero temperature effective potential, which is UV divergent. In 
the MS regularizatrion scheme, it reads 



where /U is the renormalization scale. 

The one- loop improved potential is obtained by substituting in eq. (|4,8[) the tree-level 
mass with the one-loop improved one, computed in the 4D reduced model. For g <C A, 
this simply amounts to shifting (a/R) 2 — > (a/R) 2 + XT 2 /12. We plot the thermal masses 
obtained by the one-loop naive and improved potentials in figure 19 (b) for LT ~ 1, A ~ 1 
and setting fi = 1/R, which is the scale above which our effective theory breaks down. 
The comparison between figures 19 (a) and (b) clearly indicates that the 5D model is less 
sensitive to higher-loop effects than the 4D reduced theory, although both models are not 
much sensitive to the daisy resummation procedure for a > 0.05. 

Let us now comment on the daisy diagram contributions to the Higgs effective potential 
in 5D gauge-Higgs unification models where no fundamental scalars are present. Fermion 
loops do not give rise to IR divergencies since the Matsubara modes are half integers 
and there is always an effective non- vanishing thermal mass. Moreover, chiral symmetry 
implies that fermion mass corrections are always proportional to their tree-level mass and 
thus they vanish in the IR limit of vanishing mass. There is then no need of resumming 



t f°° r 

V(T, a) = V(0, a) + — dxx 2 log 1 - e"V * 2 + ^ RT ) 2 



(4.8) 



Jo 




(4.9) 
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daisy diagrams for fermions. Gauge invariance forbids the appearance of any thermal 
mass correction to the transverse polarization of gauge bosons which are also unaffected 
by the daisy resummation procedure. The Higgs field itself, H, does not have tree-level 
self-interactions and hence does not give any higher-loop contribution to its potential in 
the daisy approximation. Finally, we are left with the longitudinal components Aq of 
the gauge fields, that acquire a non- vanishing thermal mass correction (Debye mass); the 
corresponding contributions to the Higgs potential have to be daisy-resummed. We expect 
that the latter will be of the same order of magnitude as the one we have estimated in our 
simple 5D model defined by the Lagrangian 1)4. 

As it happens in the SM, we might expect that the resummation of daisy diagrams in 
our 5D model leads to a change in the cubic term of the approximate Wilson line potential 
1)2.17(1 . However, since transverse gauge bosons are not unaffected by daisy diagrams, this 
shift would only lead to a slight change of the potential, which will still predict a first-order 
phase transition. 

As final remark, note that the absence of a tree-level potential in Wilson line based 
models permits to evade Weinberg's rough argument for the necessity of the breakdown of 
perturbation theory around the phase transition [33]: since the tree- level potential appear- 
ing in generic quantum field theories is temperature independent, perturbative corrections 
(which are instead temperature-dependent) cannot give rise to a drastic change in the 
potential (as needed for a phase transition), unless perturbation theory breaks down. 17 

5 Conclusions 

We have studied various aspects of the dynamics of Wilson line phases at finite tempera- 
ture. From a more theoretical viewpoint, we have first of all shown that the one-loop effec- 
tive potential is gauge invariant, contrary to ordinary effective potentials in non-abelian 
gauge theories in 4D, including the SM. We have also pointed out that the non-local na- 
ture of the Wilson line potential, which already ensures its finiteness at zero temperature 
at all orders in perturbation theory, is responsible for its mild dependence on quantum 
corrections at finite temperature. We have established this result by computing the one- 
loop improved thermal mass — obtained by resumming certain IR divergent higher loop 
diagrams — with the naive one in a given 5D model, and showing that the improvement 
obtained by the resummation is small. This gives strong evidence that in models in which 
the order parameter is given by a Wilson line phase (models of gauge-Higgs unification), 
perturbation theory generally holds at least up to the critical temperature Tc ~ 1/L. 

17 The possibility that perturbation theory does not necessarily break down in our 5D model is also 
suggested by the value of the critical temperature Tc ~ 1/L, which is independent on any coupling 
constant. This has to be contrasted, for instance, with the SM case, in which roughly Tc ~ \/—fi 2 /X, 
where /i 2 is the Higgs mass term. 
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Of course, it would be interesting and reassuring to establish this result by performing a 
systematic study of higher loop corrections for these 5D models. 

At a more phenomenological level, we have studied how the phase transition occurs 
in Wilson line based models, showing that this is typically of first-order. We have then 
focused our attention to a particular model in which the Wilson line phase transition is 
identified with the electroweak phase transition. In this model, unless some generalizations 
are considered, the Higgs mass is always predicted to be too low: Mjj < Mw/2. In this 
range, the transition is strongly of first order, similarly to what happens for the SM in the 
same regime of parameters. 

Introducing bulk fermions in large representations of the SU(3) W gauge group allow 
for reasonable Higgs masses. In this case, the first order phase transition is very weak and 
most likely unable to give a successful electroweak baryogenesis. Another possible option 
to increase the Higgs mass is obtained by introducing large localized gauge kinetic terms. 
Interestingly enough, the first order phase transition is now considerably stronger than in 
the previous case and it might even be sufficiently strong for electroweak baryogenesis. 

We tacitly assumed in our paper that the dynamics of Wilson line phases at finite 
temperature is decoupled from that of space-time itself, namely that the mechanism sta- 
bilizing the internal direction does not alter our analysis or takes place at temperatures 
higher than the critical temperature Tq- Since the compactification scale is of the same 
order of Tc, it would be interesting to see in some concrete model of radius stabilization 
if this assumption can actually be verified. 
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A Explicit form of the potential 



A.l Bulk fields contribution 



There are various ways of computing ea. Q2.1|l for bulk fields. For generality, we will report 
in the following the result that one obtains in d non-compact dimensions. A possible 
computation is to disentangle the T = and T / terms in ea. (|2.1|) and compute the 
T-dependent term as usually done in 4D for each KK state. In this case, one simply gets 
(see i.e. [HI US]) 



V(T,qa) = V(T = 0,qa) + (- 



+00 

k = — cx 



d d p 



(27T) 
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(A.l) 
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where p 2 denotes the d dimensional non-compact momentum square. 

Other useful and more explicit forms of the potential are obtained by using the relation 
tr \ogA = - J °° dt/t tr e~ tA in eq. (|2.1|) . valid up to irrelevant constants. Subsequently, one 
can use the Poisson resummation formula 
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1 



+00 
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(A.2) 



The above formula can be applied in ea. ()2.1|) to the KK modes k, to the Matsubara ones m, 
or to both, resulting in different but equivalent ways of computing the effective potential. 
In the three cases, one gets, respectively: 



V(T 2T(-) 2 ^ 1 ^ ~ cos(2vrfcg a ) 
V(T,qa) = m 2^ Z> — B ^ 
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where 



Bd+i(x) = x 2 Kd+i (x) 
2 2 

and K n are the modified Bessel functions. By Poisson resumming over m, we get 



V(T,qa) = (-) 2 "" 1 



L 



cos(2irkqa) 



E 



Bd+2(LMk) + 



(A.3) 



(A.4) 



(A.5) 



00 00 



1 E Eh 

k=— 00 m=l 



2rhri 



mV27T 



d+1 



Bd+i 



( M k m 
¥\ T 



where a Poisson resummation over k has been performed in the first term (m = 0). Finally, 
by Poisson resumming in both indices, one finds 
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(A.6) 



In eas. ()A.3|) . (|A.5|) and ()A.6j) an irrelevant a-independent term has always been omitted. 
All these relations also hold for vanishing bulk mass term, M = 0, by noticing that 
lim x ^o B a {x) = 2 a ~ 1 r(a), where a is any positive integer or half-integer number. 



A.2 Boundary-bulk fields contribution 

The contribution to the Higgs effective potential of bulk and boundary fields mixed as in 
ea. (|3.3j) is substantially more involved than that of purely bulk fields. One has essentially 
to compute V(T,qa) directly from ea. ()2.1j) . where M| are the mass eigenvalues of the 
bulk-boundary system. It turns out that the full contribution can be written as a sum of 
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a purely bulk contribution and a remaining boundary term. For the matter Lagrangian 
(|3.3ft . the bulk contribution is that given by the bulk fermions \& and ^ with e\ = e2 = 0, 
whereas the boundary contribution is a simple generalization of the one found in ref. [Hj 
at T = 0: 
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(A.7) 



In ea. (|A.7|l . A = ML/2 and ei^ = \/Zej/2 are dimensionless parameters defined from the 
bulk mass M and mixing terms ej appearing in ea. (|3.3[) . 



x = x + x and x A = (x + x + A ), 



(Ai 



and xo = (m + l/2)irLT is the properly rescaled Matsubara frequence for a fermion. The 
functions /o,i are defined as in ref. jSJ: 
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(A.9) 
(A.10) 



In presence of localized gauge kinetic terms, the gauge field contribution is given by 
the sum of the purely bulk term and an additional boundary term. The latter, for a field 
which couples diagonally with the Higgs field with charge q, can be written as 
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with x defined as in ea. (|A.8|) . but with xq = mirLT 
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